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Within the framework of density functional theory we investigate the nature of magnetism in 
various families of Fe-based superconductors, (i) We show that magnetization of stripe-type anti- 
ferromagnetic order always becomes stronger when As is substituted by Sb in LaOFeAs, BaFe2As2 
and LiFeAs. By calculating Pauli susceptibilities, we attribute the magnetization increase obtained 
after replacing As by Sb to the enhancement of an instability at (ji, n). This points to a strong con- 
nection between Fermi surface nesting and magnetism, which supports the theory of the itinerant 
nature of magnetism in various families of Fe-based superconductors, (ii) We find that within the 
family LaOFePn (Pn=P, As, Sb, Bi) the absence of an antiferromagnetic phase in LaOFeP and its 
presence in LaOFeAs can be attributed to the competition of instabilities in the Pauli susceptibility 
at (it, 7r) and (0,0), which further strengthens the close relation between Fermi surface nesting and 
experimentally observed magnetization, (iii) Finally, based on our relaxed structures and Pauli 
susceptibility results, we predict that LaOFeSb upon doping or application of pressure should be 
a candidate for a superconductor with the highest transition temperature among the hypothetical 
compounds LaOFeSb, LaOFeBi, ScOFeP and ScOFeAs while the parent compounds LaOFeSb and 
LaOFeBi should show at ambient pressure a stripe-type antiferromagnetic metallic state. 

PACS numbers: 74.70.-b 74.25.Ha 74.25. Jb,71.15.Mb,71.15.Pd 



I. INTRODUCTION 

After the discovery of the first high-T c iron-based su- 
perconductor LafOi^F^jFeAsi (denoted as 1111 com- 
pound) , the superconducting transition temperature was 
rapidly raised up to 55 K with substitution of La by 
Sm^. While various other families of Fe-based super- 
conductors were reported afterwards, like the 122 com- 
pounds AEFe 2 As 2 (A J B=Ca,Sr,Ba^, the 111 com- 
pounds AFeAs (A=Li,Na)i£~— and the 11 compounds 
FeCh (C/i=Se,Te)«i— , the transition temperature has 
always been lower than the highest one observed in the 
1111 systems. Besides the continuous experimental at- 
tempts to pursue higher superconducting transition tem- 
peratures in the Fe-based compounds and deeper under- 
standing of high-T c superconductivity^ 9 -, a great effort to 
understand the origin of their phase diagram has also 
been made on the theoretical side. 

Although it is widely believed that magnetically medi- 
ated rather than phonon-mediated pairing dominates the 
superconducting state due to its proximity to a stripe- 
type antiferromagnetic phased—, the origin of the mag- 
netism is still highly under debate. Some of the ex- 
perimental work and theoretical studies based on den- 
sity functional theory (DFT) support an itinerant sce- 
nario of magnetism due to the fact that the electron and 
hole sheets of the Fermi surface are nearly neste d 23 ' 38 " — 
and correlation effects are not very strong, resulting in 
a metallic state of the parent compounds^—. In con- 
trast, some authors favor a localized pictur o 21 ' 34 ' 35 ' 49 " — 
since DFT calculations fail to reproduce the experimen- 
tally observed band splitting in the stripe-type antifer- 
romagnetic phased. In these studies the observed small 
magnetic moment is attributed to highly frustrated su- 



perexchange interactions which explain the observed low 
energy spin excitations well 5 - 5 -. Apart from the oppos- 
ing viewpoints above, various other interpretations co- 
exist, such as those that propose that the magnetism 
could come from the local Hund's rule coupling 5 ^ or from 
the coexistence of localized and itinerant electrons^— . 
A recent LDA+U calculation explains the small mag- 
netic moment by formation of magnetic multipoles 62 
while LDA+DMFT (dynamical mean field theory) cal- 
culations in conjunction with angle resolved photoemis- 
sion (ARPES) experiments suggest that involvement of 
non-local fluctuations may be crucial 6 -^. Further, a recent 
DMFT calculation stressed the importance of the inter- 
play between frustrated and unfrustrated bands within a 
two-band Hubbard model at half-filling 6 - 3 . Among these 
theories, it is presently hard to decide which is the most 
promising for the observed magnetism in the parent com- 
pounds of the iron-based superconductors. 

In fact, one of the most popular theories mentioned 
above, the itinerant scenario of magnetism, was recently 
substantially challenged both by experiment^ and den- 
sity functional theory calculations 65 . On the one hand, 
an ARPES study on the 11 compound Fei +x Te shows no 
evidence of a Fermi surface nesting at (ir, 0) 6 - 4 - while mag- 
netic order with such a wave vector is detected by neutron 
scattering 6 - 6 -. On the other hand, a DFT calculation on 
LaOFeAs, BaFe2As2 and LiFeAs based on a pseudopo- 
tential method 65 reveals that the magnetic moments are 
enhanced in all compounds by replacing As with Sb while 
the Fermi surface nesting with a nesting vector of (tt, tt) 
is found with this substitution to be enhanced only in 
the 1111 compound but suppressed in the 122 and 111 
compounds. These studies question altogether the appli- 
cability of the theory of itinerant magnetism to the 11, 
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122 and 111 systems. 

The former discrepancy in the 11 compounds was soon 
resolved by a new DFT calculation based on the full 
potential linear muffin tin orbital (FPLMTO) method 
which reconciles the theory of itinerant magnetism with 
the existing experiments on Fei+ x Te~. It shows that, 
while the Fermi surface is nested at (jr, tt) in the undoped 
FeTe as in other iron-based superconductors, doping with 
0.5 electrons due to the excess of Fe in Fei +x Te leads to 
a strong (ir, 0) nesting of the Fermi surface which cor- 
responds to the observed magnetic ordering. However, 
up to now, the second question of enhanced Fermi sur- 
face nesting in 1111 versus suppression in 122 and 111 
compounds when As is substituted by Sb still remains. 

In this work, by applying Car-Parrinello molecular dy- 
namics — based on a projector augmented wave (PAW) 
basis^, we will show that, in contrast to the results from 
a pseudopotential method^, the magnetic moment and 
the Pauli susceptibility at (tt, tt) are simultaneously en- 
hanced when As is replaced by Sb in LaOFeAs, BaFe2As2 
and LiFeAs, which strongly suggests that Fermi sur- 
face nesting is closely related to the magnetic moment 
strength, and the itinerant scenario of nesting-driven 
magnetism is still valid in the 111, 122 and 1111 com- 
pounds. By further comparing the Pauli susceptibilities 
of LaOFePn with Pn=P, As, Sb, and Bi, we argue that 
the absence and the presence of magnetism at ambient 
pressure in LaOFeP and LaOFeAs, respectively, origi- 
nate from the competition between the instabilities of 
the susceptibility at (0,0) and (tt,tt), which again in- 
dicates the importance of Fermi surface nesting for the 
description of magnetism. We predict that a stripe- 
type antiferromagnetic metallic state should be present 
in the hypothetical compounds LaOFeSb and LaOFeBi. 
Finally, we study the structural and magnetic proper- 
ties of 1111 compounds including REOFeAs (RE^Ce, 
Nd, Sm), LaOFePn (Pn=P, As, Sb, Bi) and ScOFePn 
(Pn=P, As) and predict that LaOFeSb could be a su- 
perconductor with the highest transition temperature 
among these compounds. 



if structure optimizations are performed within non-spin 
polarized calculations as has been frequently pointed out 
in the literature,2M 2 . i Z2 

We used high energy cutoffs of 408 eV and 1632 eV 
for the wave functions and charge density expansion, 
respectively. The total energy was converged to less 
than 0.01 meV/atom and the cell parameters to less 
than 0.0005 A. Part of our results are double-checked 
by the full potential linearized augmented plane wave 
(FPLAPW) method as implemented in the WIEN2k 
code^i and full potential local orbital (FPLO) method 72 . 
Results are consistent among these methods. Through- 
out the paper, the Perdew-Burke-Ernzerhof generalized 
gradient approximation (GGA) to DFT has been used if 
not specified otherwise, and comparisons with the results 
from local density approximation (LDA) were also per- 
formed. In order to determine if the Fermi surface nesting 
is the driving force for the low-temperature stripe-type 
antiferromagnetic ordering, we calculate the q-dependent 
Pauli susceptibility at uj—0 and Fermi surface cuts at 
different k z planes without magnetization. These cal- 
culations were performed with the WIEN2k code using 
RKmax = 7. While 40000 k points in each k z plane are 
used in calculating Fermi surface cuts, a three dimen- 
sional grid of 128 x 128 x 128 k and q points and the 
constant matrix element approximation are employed for 
the susceptibility. For the calculations including Ce, Nd, 
and Sm atoms in the nonmagnetic phases, we apply the 
open core approximation for the localized / electrons. All 
calculations were performed in the scalar relativistic ap- 
proximation, which usually provides a good description of 
structural properties even for heavier elements^. Thus, 
spin-orbit coupling, which could be potentially relevant 
especially for Bi compound s 74 i 75 , is neglected in the cal- 
culations for the valence electrons. However, since most 
of the weight of the Bi 6p states is well below the Fermi 
energy and irrelevant to magnetic ordering, we would ex- 
pect only minor modifications for the resulting Fermi sur- 
faces and susceptibilities. 



II. METHOD 

Throughout this paper, the Car-Parrinello^ projector- 
augmented wave^S method is employed to optimize the 
lattice parameters and internal atomic positions. These 
optimized structures are then used for all subsequent 
electronic structure calculations unless stated otherwise. 
4x4x4 k-points and doubled (y/2 x \/2 x 1) unit cells 
with stripe-type antiferromagnetic order are used when 
relaxation of all lattice and electronic degrees of free- 
dom is performed. We use time steps of 0.12 fs and fric- 
tion to cool the systems to zero temperature. Note that 
the structure optimization is performed in the magnetic 
phase. As we shall show below, whenever experimental 
structures are available, our optimized structures com- 
pare well with the experimental ones. This is not the case 



III. PAULI SUSCEPTIBILITIES AND 
MAGNETISM IN 1111, 122 AND 111 
COMPOUNDS 

As is pointed out in SectionU a DFT calculation based 
on a pseudopotential method within the SIESTA codeS^ 
reveals a disconnection between magnetism and Fermi 
surface nesting in 122 and 111 compounds, (i.e., while 
magnetization is enhanced, the Pauli susceptibility at 
q = (tt, tt) which is responsible for stripe-type antifer- 
romagnetic ordering is suppressed when As is replaced 
by Sb in BaFe2As2 and LiFeAs), and therefore ques- 
tions the scenario of an itinerant nature of magnetism. 
From our spin-polarized GGA calculations for LaOFePn, 
BaFe2P?i2 and LiFePn (Pn = As, Sb), the same trends 
in magnetism are detected as observed in Ref. EH: the 
ground states are all found to be stripe-type antiferro- 
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TABLE I: Comparison between different DFT codes of the 
structures of LiFeAs optimized within GGA. ZLi=0.3385 and 
ZAs=0.7688 are obtained from our CP-PAW calculations. 





o(A) 


6(A) 


c(A) 


m(fi B ) 


dp c _As(A) 


SIESTA 65 


5.482 


5.285 


6.190 


2.54 


2.434 


VASP 85 


5.408 


5.294 


6.237 


1.5 


2.359 


WIEN21J 6 - 








1.58 


2.382 


CP-PAW 


5.422 


5.307 


6.255 


1.56 


2.385 



previous GGA calculations, such as LiFeAs calculated 
with and WIEN2k2£, there are laxge differences 

between our results and those of Refs. ,65.87 as shown 
in Table HI Furthermore, in Table HH we show the com- 
parison between experimental and optimized structural 
data for BaFe2As2, where we find that our optimized 
structure agrees with the experimental one better than 
that from Ref. [6&J. Since the electronic band structure 
close to the Fermi level is sensitive to the lattice 
structur e 20 ' 88 , the conclusion of Ref. HH based on their 
optimized structures that there is no connection between 
Fermi surface nesting and magnetism is questionable. 
Therefore, we reinvestigate the nesting property of the 
Fermi surface. 



TABLE II: Comparison between the experimental structure 
of BaFe2As2 and the optimized structures from different DFT 
codes within GGA. The magnetic moment on each Fe is also 
shown. zas=0.6495 is obtained from our CP-PAW calcula- 
tions. 



o(A) 6(A) 


c(A) 


m(fiB) 


d Fc -As(A) 


Exp. 78 5.616 5.571 


12.943 


0.87 


2.392 


SIESTA^ 5 - 5.756 5.590 


13.04 


2.78 


2.436 


CP-PAW 5.693 5.666 


13.008 


1.98 


2.396 



FIG. 1: (Color online) Comparison of normalized static q- 
dependent Pauli susceptibilities at fixed q z = n between ar- 
senide and antimonide of (a) 111 compounds, (b) 1111 com- 
pounds and (c) 122 compounds. The normalization factors 
are the susceptibilities of the corresponding arsenide systems 
for each type of compound at qo = (0, 0, n). Please note 
that the peak position is not exactly at q w = (7r,7r,7r) since 
the electron and hole Fermi surfaces are nearly nested rather 
than perfectly nested. Here, GGA is used for the DFT calcu- 
lations. 



magnetic metallic states and magnetic moments increase 
with the substitution of As by Sb in LaOFeAs, BaFe2As2 
and LiFeAs. 

However, the magnetic moments we obtained 
are 1.6 (2.3) [i B in LiFeAs (LiFeSb), 2.0 (2.5) \i B 
in BaFe 2 As 2 (BaFe 2 Sb 2 ) and 1.8 (2.2) [i B in 
LaOFeAs (LaOFcSb) 7 - 6 -, which are notably smaller than 
those obtained from the pseudopotential metho d 65 ' 87 . 
Further comparing the optimized lattice structures, we 
find that, while our results are in good agreement with 



Fig. [T] presents the comparison of normalized q- 
dependent Pauli susceptibilities at fixed q z = ir between 
arsenides and antimonides in 111, 1111 and 122 com- 
pounds. The normalization factors are the susceptibili- 
ties of the corresponding arsenide for each type of com- 
pounds at q = (0,0,7r). We fix q z = ir because of the 
fact that spins on iron are arranged antiferromagnetically 
along the z direction as observed in experiments 7 - 7 - - — . 
However, we checked that the conclusion drawn below 
will not be changed if q z = is fixed. 

We find that in all Fe-based families, the situation 
is similar. Two peaks around qo = (0, 0, n) and q,r = 
(jr, 7r, 7r) are detected in both arsenides and antimonides, 
and the peaks around (jr, 7r, it) are always stronger than 
those around (0,0, 7r), indicating that the instability to- 
wards stripe-type antiferromagnetic ordering dominates, 
which is consistent with our spin-polarized GGA calcu- 
lations. Most importantly, we find that, in contrast to 
Ref. |H, the Pauli susceptibilities are also enhanced to- 
gether with the magnetizations when As is substituted 
by Sb in LaOFeAs, LiFeAs and BaFe2As2, which demon- 
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(a) LaOFeP (b) LaOFeAs 




FIG. 2: (Color online) Fermi surface cuts for (a) LaOFeP 
and (b) LaOFeAs along different fc z -planes, where k z = n/4, 
n e {0,1,2} in units of 2tt. The cyan (light gray) curves 
are the electron Fermi surfaces around (it, it, k z ) and the red 
(dark gray) curves the hole Fermi surfaces around (0, 0, k z ). In 
order to show the nesting properties, the hole Fermi surfaces 
at (0, 0,k z ) are shown again, shifted by (ir,n). Here we use 
GGA for the DFT calculations. 



strates a connection between Fermi surface nesting and 
magnetism and consequently strongly suggests that the 
theory of Fermi-surface-nesting-driven magnetism is still 
valid. 



IV. COMPETITION OF INSTABILITIES IN 
PAULI SUSCEPTIBILITIES IN 1111 
COMPOUNDS 



TABLE III: Comparison of representative distances d,F e -pn 
in (A) where Pn=P, As, Sb and Bi in LaOFeP, LaOFeAs, 
LaOFeSb, and LaOFeBi between different GGA optimized 
structures and experimental structures, if available. 



Exp3S CP-PAW VASP 89 SIESTAS 



LaOFeP 


2.289 


2.264 


2.232 




LaOFeAs 


2.408 


2.372 


2.357 


2.446 


LaOFeSb 




2.547 


2.50 


2.660 


LaOFeBi 




2.639 








FIG. 3: (Color online) Static q-dependent Pauli suscepti- 
bilities at fixed q z = n for (a) LaOFeP and (b) LaOFeAs. 
The corresponding values of the Pauli susceptibilities at 
qo = (0, 0, 7r) in LaOFeP and LaOFeAs, respectively, are sub- 
tracted. On top, two-dimensional contour maps are shown. 
Here, GGA is used for the DFT calculations. 



TABLE IV: Comparison of representative distances dLa-o in 
(A) in LaOFeP, LaOFeAs, LaOFeSb, and LaOFeBi between 
different GGA optimized structures and experimental struc- 
tures, if available. 



Exp. CP-PAW (ours) VASP (ref.) SIESTA 



LaOFeP 2.350 


2.344 


2.349 




LaOFeAs 2.363 


2.356 


2.369 


2.374 


LaOFeSb - 


2.375 


2.394 


2.398 


LaOFeBi - 


2.382 







In what follows we concentrate on the 1111 compounds 
and perform a comparative study among LaOFePn 
(Pn=P, As, Sb, Bi). In Table HH and [TV] we first present 
the comparisons of two representative atomic distances 
among different structures optimized within GGA and 
experimental structures, if available. Our results agree 



well with the experimental ones. 

Fig. [5] shows the calculated Fermi surface cuts for 
LaOFeP and LaOFeAs on different fc^-planes based on 
the experimental lattice structures. The figure is almost 
unchanged if we consider the optimized lattice struc- 
ture. The hole Fermi surfaces around (0, 0, k z ) are shifted 
by (-7T, tt, 0) to show the nesting properties. Shifting of 
(tt, 7t, tt) was also investigated, and we find that the nest- 
ing properties are nearly unchanged. From the figure, 
it is apparent that the Fermi surface nesting is even 
more perfect in LaOFeP than in LaOFeAs, indicating a 
stronger tendency to stripe-type antiferromagnetic order- 
ing in LaOFeP compared to LaOFeAs. However, experi- 
mentally, while a small magnetization is observed in un- 
doped LaOFeAs^l, superconductivity rather than mag- 
netic order is detected in undoped LaOFeP—. These 
observations would indicate that Fermi surface nesting 
might not be connected to magnetization. 

In order to quantify the Fermi surface nesting, we show 
in Fig. [3] the q-dependent Pauli susceptibilities at fixed 
q z = tt for LaOFeP and LaOFeAs with subtraction of the 
corresponding values at qo = (0,0, tt). While a peak in 
LaOFeP appears right at q^ = (tt, tt, tt) indicating almost 
perfect nesting properties of the Fermi surfaces, peaks 
are situated close to q,r = (tt, tt, tt) in LaOFeAs sug- 
gesting nearly nested Fermi surfaces, which is consistent 
with the Fermi surface cuts shown in Fig. [5] The most 
interesting finding in Fig. [3] is that the relative values 
X(7r, tt, tt)~x{0, 0, tt) increase from LaOFeP to LaOFeAs 
irrespective of whether the Fermi surface nesting is per- 
fect or not. While the peak at (tt, tt, tt) favors stripe- 
type antiferromagnetic ordering, the one at (0,0, tt) rep- 
resents a possible instability towards checkerboard-type 
antiferromagnetic ordering or A-type antiferromagnetic 
ordering where ferromagnetic layers are stacked antifer- 
romagnetically. The heights of these two peaks become 
closer in LaOFeP than in LaOFeAs, implying that com- 
petition between the above-mentioned two types of anti- 
ferromagnetic states becomes stronger in LaOFeP if ther- 
mal or quantum fluctuations are taken into account. Also 
spin-fluctuation mediated pairing of the superconducting 
stat u 23 ' 26 ! 30 coming from inter-band scattering around 
q?r = (f, A") tt) takes part in the competition. Eventually, 
as the two types of antiferromagnetism strongly compete 
with each other, the additional superconducting state or- 
der emerges and opens a gap, removing the high instabil- 
ity at the Fermi level and lowering the total energy. This 
could be the scenario to explain why undoped LaOFeP 
is always nonmagnetic but shows superconductivity be- 
low 3.2 K at ambient pressure. As xi^i ^ 7r ) — x(0; 0; 7r ) 
increases beyond a critical value, the stripe-type antifer- 
romagnetic ordering prevails over the checkerboard-type 
one and the pairing. This scenario may apply to the 
low-temperature magnetic phase of LaOFeAs. Further- 
more, we have calculated the total energies of checker- 
board and stripe-type antiferromagnetic phases for both 
LaOFeP and LaOFeAs. We found that the stripe-type 
antiferromagnetic phases are the ground state in both 




FIG. 4: (Color online) Static q-dependent Pauli suscepti- 
bilities at fixed q z = n for (a) LaOFeSb and (b) LaOFeBi. 
The corresponding values of the Pauli susceptibilities at 
qo = (0, 0, tt) in LaOFeSb and LaOFeBi, respectively, are sub- 
tracted. On top, two-dimensional contour maps are shown. 
In the DFT calculations GGA is used. 



cases and the energy difference between the two phases 
is smaller in LaOFeP than in LaOFeAs, which is consis- 
tent with the trend of x(7r, tt, tt)~x(0, 0, tt). 

In Fig. HJ we display the q-dependent Pauli suscepti- 
bilities at fixed q z = tt for the hypothetical compounds 
LaOFeSb and LaOFeBi. The corresponding values of 
the Pauli susceptibilities at qo = (0, 0, tt) in LaOFeSb 
and LaOFeBi, respectively, are again subtracted. Note 
that x(7r,7r, tt) — x(0,0, tt) in LaOFeSb and LaOFeBi is 
even larger than that in LaOFeAs (see Fig. [3]), indicat- 
ing that the instability towards stripe-type antiferromag- 
netic ordering could win the competition between dif- 
ferent instabilities in these two compounds. It is also 
interesting to note that the peak at qo = (0, 0, tt) be- 
comes flatter when we go from LaOFeP to LaOFeBi. 
While the flatness of the peak can be associated with 
a larger number of different magnetic structures lying 
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FIG. 5: (Color online) Static Pauli susceptibilities at qo = 
(0, 0, 71") and = (tt, tt, tt) for LaOFePn with Pn=P, As, Sb, 
Bi. Here we use GGA for the DFT calculations. 

within a small energy window, the stronger peaks around 
qjr = (7r,7r, 7r) compared to qo = (0,0, it) in all the 1111 
compounds we studied makes other magnetic orderings 
besides the stripe-type antiferromagnetic one less favor- 
able. Combining the results from spin-polarized GGA 
(LSDA) calculations where magnetic moments on each 
iron are given as 2.2 (1.3) and 2.4 (1.8) hb in LaOFeSb 
and LaOFeBi, respectively, we predict that the ground 
states of these two compounds at ambient pressure with- 
out doping should show stripe-type antiferromagnetic or- 
der although spin-polarized GGA (LSDA) calculations 
overestimate the magnetic moments. 

The increase of the magnetic moment from As to Sb to 
Bi can be understood from Fig. [5] where Pauli suscepti- 
bilities at qo = (0, 0, 7r) and q^ = (tt, tt, it) for LaOFePn 
with Pn—P, As, Sb, Bi are explicitly shown. While the 
increasing absolute value of XPnl 71 ", tt, tt) as Pn changes 
from P to Bi is probably responsible for the increasing 
magnetic moment, the difference between xp n (TT,TT,Tr) 
and XPn(0,0,7r) dominates the possible competition be- 
tween different ordered states, which again implies a 
strong relation between Fermi surface nesting and mag- 
netism. Due to the fact that LaOFeAs is not a super- 
conductor without doping or application of pressure, we 
argue that possible superconducting states in LaOFeSb 
and LaOFeBi can only occur under doping or application 
of pressure. 

In Fig. © we display the DOS for LaOFeAs, LaOFeSb 
and LaOFeBi calculated within spin-polarized GGA cal- 
culations. It is shown that in all three cases the DOS at 
the Fermi level remains finite, suggesting that undoped 
LaOFeSb and LaOFeBi are stripe-type antiferromagnetic 
metals at ambient pressure as is the case of LaOFeAs 
or other iron-based superconductors. Our results from 
spin-polarized GGA calculations cannot corroborate the 
arguments introduced in Ref. [9l| where it is argued that 
LaOFeSb and LaOFeBi could be antiferromagnetic Mott 
insulators, though the tendency towards larger Fe mag- 




-1 1 

energy (eV) 

FIG. 6: (Color online) Total density of states (DOS) for 
LaOFeAs, LaOFeSb and LaOFeBi in the low-temperature or- 
thorhombic phase with stripe-type antiferromagnetic order 
calculated by spin-polarized GGA. Also shown for compari- 
son is the DOS for nonmagnetic LaOFeP calculated by GGA. 
Relatively high DOS at the Fermi level is present in LaOFeP, 
indicating the possible instability with respect to supercon- 
ductivity in the absence of magnetic order. 



netic moments from the As to Bi compounds is observed. 
Such a conclusion is also confirmed by our LSDA calcula- 
tions. In Fig. [SI the DOS for LaOFeP is also shown. Since 
LaOFeP is nonmagnetic at low temperature, only the 
non-spinpolarized GGA result is presented. It is found 
that the DOS at the Fermi level remains relatively high, 
which indicates a possible instability with respect to su- 
perconductivity in the absence of magnetic order with 
the opening of a superconducting gap and the lifting of 
the degeneracy at the Fermi level. 



V. NEW SUPERCONDUCTOR CANDIDATES 

Experimentally, the highest recorded superconducting 
transition temperature has been observed in the 1111 
compounds. It is therefore tempting to find ways to pre- 
dict T c for hypothetical 1111 compounds. One promis- 
ing route is to consider a phenomenological relation be- 
tween T c and the x(7r, tt, tt) — x(0, 0, tt) of the parent com- 
pound. As we know, with doping or application of pres- 
sure, the x( n > 7r ) instability is suppressed, resulting in 
the disappearance of stripe-type antiferromagnetic order. 
However, strong inter-band scattering with a wave vec- 
tor around (tt, tt, tt) compared to intra-band scattering 
with a wave vector around (0,0, tt) remains, which leads 
to a superconducting state. Therefore we argue that the 
larger the relative value of x(7r, tt, 7r)-x(0, 0, tt) in the par- 
ent compound is, the stronger the antiferromagnetic spin 
fluctuation after suppression of magnetic order, and thus 
the higher the superconducting transition temperature 
will be. 
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FIG. 7: (Color online) Prediction of superconducting tran- 
sition temperatures T c from a phenomenological relation be- 
tween T c and x( 7r i ""j 7r ) ~ x(0j 0, 7r) for the parent compounds 
of the hypothetical 1111 compounds LaOFeSb, LaOFeBi, 
ScOFeP and ScOFeAs. The phenomenological relation was 
determined by calculating x(n, it, it) — x(0, 0, it) for several 
1111 compounds LaOFeP, LaOFeAs, CeOFeAs, NdOFeAs, 
and SmOFeAs by GGA where T c 's and lattice structures are 
given experimentally. The x( n i 17 1 7r ) — X(0, 0; 7r ) f° r LaOFeP 
and LaOFeAs calculated by GGA from DFT optimized struc- 
tures are also shown for comparison. It is found that the re- 
sulting x( 7r i 1") 7r ) — x(0, 0, 7r) for LaOFeP and LaOFeAs based 
on our optimized structures is only slightly underestimated 
compared to that calculated from experimental structures, 
which shows that the result depends only weakly on our struc- 
ture optimization. 
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FIG. 8: (Color online) Prediction of superconducting tran- 
sition temperatures T c from a phenomenological relation 
between T c and the ratio of atomic distances dFe-Fc and 
dFc-As for the hypothetical 1111 compounds LaOFeSb, 
LaOFeBi, ScOFeP and ScOFeAs. A constant ratio of 
ro = do,Fc-Fc/do,Fc-As for distances do,Fe-Fc, do,Fc-As in a per- 
fect FeAs4 tetrahedron is subtracted. The phenomenologi- 
cal relation is determined by taking into account the same 
compounds LaOFeP, LaOFeAs, CeOFeAs, NdOFeAs, and 
SmOFeAs as in Fig. where T c and lattice structures are 
given experimentally and doing a linear fit as in Ref. |95| . 
The |dFc-Fo/dFe-A s - r | for LaOFeP, LaOFeAs, LaOFeSb, 
and LaOFeBi based on optimized structures within GGA and 
LDA are also shown for comparisons. Comparing the results 
for LaOFeP and LaOFeAs from GGA and LDA optimiza- 
tions with those from experiments, we find that the optimized 
structures from GGA are more consistent with the experimen- 
tal one. 



In Fig. [7J we plot xi^^^n) — x(0,0, it) versus T c for 
the hypothetical 1111 compounds LaOFeSb, LaOFeBi, 
ScOFeP and ScOFeAs. The phenomenological relation 
between T c and x(7r,7r, it) — x(0,0,7r) is determined by 
first calculating x( 7r , 7r 7 7r ) — x(0>0, 7r) for several typ- 
ical 1111 compounds LaOFeP^ 2 ., LaOFeAsiH, Ce- 
OFeAsH^, NdOFeAaSkSI, and SmOFeAs^^, where T c 
and lattice structures are given experimentally, and then 
fitting the data by an exponential growth function of 
X(7r,7r,7r)-x(0,0,7r) = 2.86 + 0.28*exp(T c /18.14). With 
this relation, Tc for the compounds which have not yet 
been experimentally reported is predicted by optimizing 
the lattice structure and calculating x(tt, tt, it) — x(0, 0, n) 
from DFT calculations. From Fig.[7j we find that the re- 
sulting x( 7r , 7r , 7r ) — x{0,0,^) for LaOFeP and LaOFeAs 
based on our optimized structures is only slightly under- 
estimated compared to that calculated from experimen- 
tal structures, which shows that the result depends only 
weakly on our structure optimization. Among the four 
1111 compounds, this procedure shows that LaOFeSb 
can give the highest T c around 60 K which is above the 
highest recorded T c of 55 K in SmOFeAs. 

An alternative procedure to predict T c phcnomcnolog- 
ically is based on the fact that the physical properties 



of Fe-based superconductors strongly depend on the As 
position. Tiny shiftings of the As position away from 
or closer to the iron plane will significantly change the 
band structure around the Fermi leve l 20 ' 88 . Therefore, 
similar to Refs. HH and [lOOl we plot in Fig. U T c versus 
the absolute value of the ratio between atomic distances 
d F c-Fo and d Fc _ As while subtracting r = c?o,Fc-Fc/do,Fc-As 
where do,Fc-Fc and do,Fe-As denote the distances in a 
perfect tetrahedron formed by four nearest neighbor As 
atoms surrounding one Fe atom. Similar to the first 
scheme, we determine the phenomenological relation be- 
tween IdFe-Fe/^Fe-As — tq\ and T c by taking into account 
the same compounds LaOFeP, LaOFeAs, CeOFeAs, Nd- 
OFeAs, and SmOFeAs as in the first scheme where T c 
and lattice structures are given experimentally and do- 
ing a linear fit as in Ref. 95. With this relation, the T c for 
LaOFeSb, LaOFeBi, ScOFeP and ScOFeAs is predicted 
based on the optimized structure. 

Fig. M presents the results based on the structures op- 
timized within both GGA and LDA. Comparing the re- 
sults for LaOFeP and LaOFeAs from GGA and LDA op- 
timizations with those from experiments, we find that 
the optimized structures from GGA are more consis- 



8 



tent with the experimental one. According to the re- 
lation we fitted, LaOFcSb always gives the highest T c of 
57.5 (50.6) K for LDA (GGA) optimizations, respectively, 
among the four 1111 compounds we studied. Combining 
the two presented phenomenological prediction schemes, 
we clearly obtain that LaOFeSb would be upon doping 
or under pressure a good candidate for superconductivity 
with highest T c and it would be very interesting to see it 
synthesized. 

VI. CONCLUSIONS 

In the present work we studied the physical proper- 
ties of LaOFePn, BaFe2-Pn2 and LiFePn with Pn=As 
and Sb. Our results support the validity of the itin- 
erant nature of magnetism in these compounds where 
magnetization is closely related to Fermi surface nest- 
ing. Furthermore, we concentrated on the 1111 com- 
pounds LaOFePn with Pn~P, As, Sb, and Bi. We 
found that the increase of the magnetic moment in the 
undoped compounds with Pn varying from P to Bi is 
due to the increasing instability of the Pauli suscep- 
tibility at q^ = (7T,7r,7r) and the decreasing competi- 
tion to the instability at q = (0,0, n). The supercon- 
ducting state appearing in undoped LaOFcP at ambi- 
ent pressure is ascribed to the strong competition be- 
tween the instability at q w = (tt, tt, tt) and qo = (0, 0, tt). 
Thus, together with the investigation of the DOS in the 



low temperature phase, we argue that the hypotheti- 
cal compounds LaOFcSb and LaOFcBi are antiferromag- 
netic metals at ambient pressure without doping. The 
results for LaOFePn again strongly imply that Fermi 
surface nesting plays a dominating role in the physical 
properties of the 1111 compounds. Finally we consider 
two phenomenological relations to predict the supercon- 
ducting transition temperature T c for the hypothetical 
1111 compounds and predict that LaOFeSb would be a 
possible candidate for a superconductor with a higher 
T c than presently recorded for the known Fe-based su- 
perconductors. Combining the fact that Fermi surface 
nesting dominates the physics in the 122 compounds, 
1111 compounds, 11 compounds, and 111 compounds, 
we argue that magnetism in iron-based superconductors 
is strongly influenced by their itinerant nature. However, 
from our study, the localized scenario is not ruled out and 
may also play an important role in the physics of iron- 
based superconductors. Furthermore, while in our study 
we emphasize the role of the states at the Fermi level and 
accordingly the nesting property of the Fermi surface on 
the itinerant nature of magnetism, the significant contri- 
butions to the finite moment of itinerant magnetism from 
the states in the vicinity of the Fermi level should not be 
ignored. 
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